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Abstract 

There is an approach due to Bazhanov and Reshetikhin for solving integrable RSOS 
models which consists of solving the functional relations which result from the trun¬ 
cation of the fusion hierarchy. We demonstrate that this is also an effective means of 
solving integrable vertex models. Indeed, we use this method to recover the known 
Bethe Ansatz solutions of both the closed and open XXZ quantum spin chains with 
U(l) symmetry. Moreover, since this method does not rely on the existence of a pseu- 
dovacumn state, we also use this method to solve a special case of the open XXZ chain 
with nondiagonal boundary terms. 


1 Introduction 


There are several well-known methods of deriving the Bethe Ansatz (BA) solution of inte¬ 
grate vertex models: the coordinate BA (T|, |2|. ||, the T — Q approach ||, the algebraic BA 
||, |[. | 6 ]], the analytic BA [[7|, and the functional BA ||. We present here yet another method, 
which entails solving the functional relations which result from the truncation of a model’s 
fusion hierarchy. This approach was (to our knowledge) first developed for RSOS models 
by Bazhanov and Reshetikhin in [|TD[, but until now has not been applied to vertex-type 
models. An important feature of this method is that, unlike some of the other approaches, 
it does not rely on the existence of a pseudovacuum (reference) state. 


Our primary motivation for this work comes from the long outstanding problem of solving 
the open spin-1 XXZ quantum spin chain with nondiagonal boundary terms, defined by the 


Hamiltonian 11 


N -1 


ft = 2 { ( a n a n+ 1 + <° V n+l + COsh V «+l) 


n= 1 


/ 2k 

+ sinh 77 ( coth £_cr? -|--- af — coth£ + cr 

V sinh F- 


sinli 


2k i 


sinh £ + 


-a 


N 


( 1 . 1 ) 


where a x ,a y , a z are the standard Pauli matrices, q is the bulk anisotropy parameter, , k± 
are arbitrary boundary parameters, and N is the number of spins. This model is integrable. 
Indeed, the Hamiltonian is obtained from the commuting transfer matrix [0] constructed 
with the nondiagonal K matrix found in [II], |T2[ together with the standard XXZ R matrix. 


Solving this problem (e.g., determining the energy eigenvalues in terms of roots of a 
system of Bethe Ansatz equations) is a crucial step in formulating the thermodynamics of 
the spin chain and of the boundary sine-Gordon model Jl2|. Moreover, this problem has 
important applications in condensed matter physics and statistical mechanics. 


A fundamental difficulty is that, in contrast to the special case of diagonal boundary 
terms (i.e., k± = 0, in which case H has a U( 1) symmetry) considered in [|], |], a simple 
pseudovacuum state does not exist (e.g., the state with all spins up is not an eigenstate of 
the Hamiltonian). 


Some progress on this problem was made recently in |13|. |1 4}| . Namely, for bulk anisotropy 
value r] = , p — 1 , 2 ,..., (and hence q = e ,? is a root of unity, satisfying q p+1 = — 1 ), 

an exact (p + l)-order functional relation for the fundamental transfer matrix was proposed. 
The key observation is that the fused spin-^yp transfer matrix can be expressed in terms 
of a lower-spin transfer matrix, resulting in the truncation of the fusion hierarchy, [j The 
1 This is distinct from the observation due to Belavin et al. pl| nR that, for the special case of quantum- 


1 

















simplest case p — 1, which corresponds to the XX chain, is analyzed in [13 . 

Although sets of equations for the eigenvalues of the transfer matrix were found in [[14] 


these equations do not have the standard Bethe Ansatz form, and they become increasingly 
complicated as the value of p increases. Moreover, one would like to solve the model for 
general values of ??, i.e., not just for the discrete values corresponding to roots of unity. 

We also achieve here some progress on these questions. In particular, from the func¬ 
tional relations in [|14] , we obtain standard Bethe Ansatz equations for the transfer matrix 
eigenvalues for general values of p, albeit only for the special case 


k+ = k- = k ± 0 , £+ = i- = £ , 


N = odd. 


( 1 . 2 ) 


Unfortunately, we have not yet succeeded to obtain corresponding results for general values 
of the boundary parameters. 

The outline of this article is as follows. In Section |2|, we consider as a warm-up the case of 
the closed XXZ chain, and provide a new derivation of the well-known Bethe Ansatz solution. 
In Section [3]. we turn to the open XXZ chain. Using the functional relations proposed in 
14 1, we first recover the solution of Alcaraz et al. [|| and Sklyanin [f| for the diagonal case 
k± = 0, and we then give the solution for the nondiagonal case (|1.2|) . We conclude with a 
brief discussion of our results in Section [|. 


2 The closed chain 


The closed (i.e., with periodic boundary conditions) spin-4 XXZ quantum spin chain is 
defined by the Hamiltonian 

1 N 

H =~Y1 ( a n a n+ 1 + a n a n+ 1 + COsh V «+l) , (2-1) 

n =1 


where cqv+i = (7\■ As noted in the Introduction, there are various methods of deriving the 
Bethe Ansatz solution of this model. As a warm-up for the open-chain problem, we now give 
another derivation of this solution, which involves solving the model’s functional relations. 


In Sec. |2.1| we derive the functional relations, and then in Sec. |2.2| we proceed to solve them. 


group symmetry (i.e., k± = 0, £± —> oo), the fused transfer matrix t^\u) vanishes after quantum group 
reduction. 
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2.1 Functional relations 


In this subsection, we begin by briefly reviewing the construction of the (fused) transfer 
matrices of the closed XXZ chain. We then recall the so-called fusion hierarchy which these 
transfer matrices obey. Finally, we give an identity which truncates the fusion hierarchy and 
leads to the desired functional relations. 

The fundamental spin-(|, |) XXZ R matrix is given by the 4x4 matrix 

/ sinh(w + rj) 0 0 0 

0 sinh u sinli rj 0 

R(u) = 

0 sinh rj sinh u 0 

\ 0 0 0 sinh(u + rj) 

where rj is the anisotropy parameter. It is a solution of the Yang-Baxter equation 


( 2 . 2 ) 


R 12 (u - v ) R 13 (u) R 23 (v ) = R 23 (v) R 13 (u) R 12 (u - v) . (2.3) 

(See, e.g., 0, ITsfl.) The fused spin-(j, |) R matrix (j — \ , 1, |,...) is given by |l^] 
R(i...2j)2j+i(u) = Pi 2j Ri :2j+1 (u)R 2 , 2j+1 (u + rj)... R 2jt2j+1 (u + (2 j - 1 )rj)Pl 2j . (2.4) 

The (undeformed) projectors are defined by 

pt m =-,Y. , ~ ±i y v -’ < 2 ' 5 ) 

ml z —^ 


(7 

where the sum is over all permutations a = (fly ,... cr m ) of (1,... , m), and V a is the permu¬ 
tation operator in the space ®1 =l C 2 . For instance, 


P\2 ~ 2 ^ + ^ 12 )’ 

1*123 = “(I + V 23 P\ 2 + Pl 2 P 2 3 + Fl 2 + V 2 3 + F 13 ) , ( 2 . 6 ) 

6 

where I is the identity matrix. 

The closed-chain transfer matrix t^\u), which is constructed using a spin-j auxiliary 
space, is defined by 


= tr, ..2, / 1...2/. ( 11 ), 


where the fused monodromy matrix is defined by 

= R i J,. Vi " J - - ■ , 


(2.7) 


( 2 . 8 ) 
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and N corresponds to the number of spins of the chain. One can show that 


= Pi ;_ 2 i Ti(u)T 2 (u + rj) ■ ■ ■ T 2j {u + (2 j - 1 )v)Pl. 2j ■ 

These transfer matrices constitute commutative families 

[t^\u ), l k \v)] = 0 . 

The fundamental transfer matrix t{u) = t^\u) contains the Hamiltonian (|2.1|), 


d 

H (x — logt(u) 


u =0 


+ const ., 


and has the periodicity property 


t{u + in) = (—1)^ t(u) . 


The fusion hierarchy for the XXZ chain is given by [17], |21 


t°\u) t^\u + 2jrj) = 5(u+ (2 j - 1 )rj) t u ?\u) + t (j+ i\u ), j = - , 1, - , 
where t^°\u) = I, and the quantum determinant [17], |2(J 5{u) is given by 

N 


where 


5(w) = tris P 12 Ti(u)T 2 (u + rj) = (-C(u + 7])) , 


C(u) — — sinh(u + rj) sinh(w — rj) . 


The key fact in deriving the functional relations is that for anisotropy values 

in 


V = 


p+ 1 


P = 1,2,. 


the fused transfer matrices satisfy the identity 

l^lu) = (—1 ) N S(u — rj) t {E ^\u + rj) + (1 + (—1 ) N )v{u) N F 

where 


n(u) = 


C(n) 


sinh(u + krj) 


i ay 

- ) sinh((p + l)w) , 


fc =0 


CM \2 


and 


N 


F = U.° 

k=1 


k * 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

, (2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 
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(2.19) 


















The remarkable result (2.17), to which we refer as the “truncation identity,” follows directly 


from Eq. (4.13) in Ref. |lj], which relies on the quantum-group construction |22| of higher- 
spin R matrices. 

The fact that the spin- 241 transfer matrix can be expressed in terms of a lower-spin 
transfer matrix leads to the truncation of the fusion hierarchy, which in turn leads to a 
(p + l)-order functional relation for the fundamental transfer matrix. For instance, for the 
case p = 2, Eqs. (|2.13|) and ( |2.17| ) lead to the third-order functional relation 

t(u)t(u + rf)t{u + 2 rf) — S(u)t(u + 2 rf) — S(u + r))t(u) 

— (—1 ) N S(u — rj)t(u + rf — (1 + (—1 ) N )5(u — if)v(u) N F = 0 . (2.20) 

Similar higher-order functional relations have been obtained for RSOS models |2|, ^3], [10| and 
for the 8-vertex model |24|. We emphasize that, contrary to the commonly-held misconcep¬ 
tion (see, e.g. |23|), the fusion hierarchies of vertex models do truncate, for the rj values 

(B 

The commutativity relation ( 2.10|) with j = k = | and the fact [F ,t(u)] = 0 imply that 
t(u) and F can be simultaneously diagonalized, 

t(u) |A (±1) ) = A (±1) (m)|A (±1) ) , 

F|A (±1) ) = ±|A (=tl) ), (2.21) 

where the eigenstates jA^ 1 ^) are independent of u. Acting on these eigenstates with the 
functional relations, one obtains the corresponding relations for the eigenvalues. 


2.2 Bethe Ansatz solution 

We now proceed to solve the functional relations for the eigenvalues of the fundamental 


transfer matrix. Following Ref. |[T0|, we observe that the functional relations for p > 2 can 
be represented in a compact form as the determinant of a (p + 1) x (p + 1) matrix: 


/ A, 


det 


(F) 


-h 

0 


where 


-h -1 
Af } 

-ho A 


0 0 
\ -Fh p _ i 0 


0 

-ho 

(F) 


0 

0 


0 

0 

-hi 

0 

0 


0 

0 

0 


0 

0 

0 


-Fh 0 

0 

0 


\ 


(F) 

P- 

~ hr) A 


-/ip-1 Ap_) 


(F) 


= 0 , 


( 2 . 22 ) 


h(u) = sinh A (u + rf ), 


(2.23) 
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hk = h(u + r)k), A[ F) = (u + rjk) , and F — ±1 now denotes the eigenvalue of the operator 

in (|2.19|) . Let (Q 0 ,Q i,... , Q p ) be the null vector of the matrix in (ED; i-e-, 


A^'Qo - h-iQi - Fh 0 Q r = 
hkQk —1 ”1“ Qk lQ/c+1 

—Fhp-iQo — hpQp- 1 + A^Qp = 
We make the Ansatz Q*, = Q(u + rjk), where 

M 

Q(u ) = JJsinh(n — 
j=i 


k= 1, ■ • • ,P-1, 




'j/; 


(2.24) 


(2.25) 


for some integer M. Eqs. ( |2.24|) imply (using Qk = (—l) M Qk+ P +i) that the eigenvalues are 
given by 


A m {u) = h{u) Q^pI!l + h{u _ rl) Q^+^ 


and 


Q(u) 


F=(~ 1 ) 


Q(u) 


M 


(2.26) 


(2.27) 


We verify that the result 


is consistent with the periodicity condition (2.12). The 


requirement that A F)[u) be analytic at u = Uj yields the Bethe Ansatz equations 

h(uj) _ Q(uj + 7]) 


h(uj - rj) Q(uj - rj) ’ 


3 = !,••• ,M. 


(2.28) 


We recognize Eqs. ( |2.23| ), (|2.25|) , ( |2.26|) and ( |2.28|) as the familiar Bethe Ansatz result for 
the eigenvalues of the transfer matrix of the closed XXZ chain. Although we have assumed 
that r) has the values ( |2.1G| ), these results are known to be true for general values of 77 . Note 
also that the approach that we have followed here does not explicitly rely on the existence 
of a pseudovacuum state. 

A more thorough analysis would also include the diagonalization of S z = ^ J2k =1 a k 
along with t(u) and F. By considering the asymptotic behavior of t(u) for u —* 00 , one can 
establish that the value of M is related to the S z eigenvalue, namely, M = A- — S z . Since 
these matters are already well-understood, and since the open-chain problem (| 1 . 1 |) lacks this 
additional U( 1) symmetry, we do not pursue this issue further. 

Finally, we point out that the result ( 2.27|) , which perhaps is less familiar, can nevertheless 
be readily obtained within the algebraic Bethe Ansatz approach ||. Indeed, as is well known, 
the eigenstates of the transfer matrix are given by 


\ui ,... ,u M ) = B(ui). ..B(u M )\ty , 


(2.29) 
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where B(u) is a certain creation-like operator, and |fl) is the pseudovacuum state with 
all spins up. It is easy to show that {F ,B(u)} = 0 and F\D) = |f2). Hence, the state 
|Mi,... , um) has the eigenvalue F = (—1) M . 


3 The open chain 


We turn now to the open chain (|1.1|), which is our main concern. Our strategy is to try to 


generalize the analysis of the preceding section. Hence, in Sec. |3.1| we review the functional 
relations, and then in Sec. |3.2| we attempt to solve them. 


3.1 Functional relations 


The fundamental spin-7 XXZ K matrix is given by the 2x2 matrix il 

K-(u) = I ' '■ '~T.I . (3.1) 


. . / sinh(£_ + m) k,- sinh 2u 

( “ ) = K- sinh 2u sinh«. - „) 


which evidently depends on two boundary parameters , k_. It is a solution of the boundary 
Yang-Baxter equation [2d| 


Rv 2 {u - v) K 1 (m) R 2 i(u + v) K 2 (v) = K 2 (v) R\ 2 [u + v) K l ( u) R 2 \(u - v ). (3.2) 

The fundamental open-chain transfer matrix t(u) is constructed, following Sklyanin’s 
recipe [Q, from the matrix R(u) ( | 2 . 2 | ), the matrix K~(u ) O. and the matrix K + {u ) which 
is equal to K~(—u — rj) with (£_,/c_) replaced by (£ + ,/c + ). Q The Hamiltonian O is 
related to the hrst derivative of the transfer matrix, 


H = 


dt{u) 


4 sinh sinh sinh 2 1 r] cosh r] du 
The transfer matrix has the periodicity property 

t(u + m) = t(u ), 

as well as crossing symmetry 

t(—u — rj) = t(u ), 

and the asymptotic behavior (for k± ^ 0 ) 

e «(2V+4)+r ) (Af+2) 


u =0 


sinh 2 77 + N cosh 2 77 
2 cosh 7 ] 


I. 


t(u) 


-K-K+- 


22N+1 


-I+... 


for 


u —> 00 . 


(3.3) 


(3.4) 


(3.5) 


(3.6) 


2 Further details about the construction of this transfer matrix can be found in 
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Functional relations for the open XXZ chain (1.1) have been proposed in Ref. [ 141. These 
relations, which follow from the fusion hierarchy |27| |28| together with the truncation identity 
for the 77 values (|2.16|), are given by 


A(u)A(u + rj)... A (77 + pp ) 

— 5(u — p)A(u + p)A(u + 2rj)... A {u + (p — 1 )p) 

— S(u)A(u + 2p)A(u + 3?/)... A (u + pr 7 ) 

— S(u + ? 7 )A(w)A(w + 3p)A(u + Arf)... A(u + pij) 

— 5{u + 2p)A(u)A(u + p)A(u + 4 rj)... A(u + pp) — ... 

— 8(u + (p — 1 )? 7 )A(w)A(-u + 77 )... A(u + {p — 2)p) 

+ ■ ■ ■ = f(u ), (3.7) 


where A (u) is the eigenvalue of the fundamental open-chain transfer matrix t(u). Further¬ 
more, the function 8(u) is now defined by 




A(«) 

C(2u + 2p) ’ 


(3.8) 


where the quantum determinant A (u) is given by 


A (77) = — [sinh(u + p + £_) sinh(u + p — £_) + n 2 _ sinh 2 (2u + 2 p)\ 

x [sinh(w + 77 + £ + ) sinh(u + p — £+) + n 2 + sinh 2 (2w + 277)] 
x sinh 2 u sinh( 2 u + 477) ((u + p) 2N , ( 3 . 9 ) 


and C( u) is defined in Eq. ( |2.15| ). Moreover, the function f(u) is given by Q 


/(«) = 


. j_)p(-W+i) 


sinh 2iV ((p + l)u) 


cosh 2 ((p + l )u + y 6 ) 


22 P (n+i) cosh 2 ((p + 1 ) 77 ) 

x \ n(u ; , «_) 77(77; -£+ , n + ) + 77(77 ; , «_) n(u ; £+ , hi + ) 


+2(—l) 7V (—7 c_tc + ) p+1 sinh 2 (2(p + 1 ) 77 )} , 


(3.10) 


3 In terms of the three functions /o, / 1 , fz used in Q, the functions S(u) and f(u) are given by 


<5(m) = 


fiju + p) 

fo(u) 


f(u) = 


fs(u) 
fo(u) ' 



















where e = 2 frac(p/ 2 ) equals 0 if p is even, and equals 1 if p is odd; and the function n{u ; £ , n) 
is defined by 

m 

n(u-,£,n) = sinh ((p + 1)(£ + u))+ ^ c P)[ K 2 ; sinh((p + l)u+(p+l — 2 1)£) , (3.11) 

i=i 

with 

k=0 

For instance, for the case p = 3, the functional relation is given by [] 

A(u)A(u + ij)A(u + 2i])A(u + 3?/) — 5(u — t/)A {u + ij)A (u + 2 rj) 

—8{u)A{u + 27/) A {u + 37 /) — 5{u + rj)A(u)A(u + 3t/) — 5{u + 2rj)A{u)A{u + 77 ) 

+< 5 ( 77 ) 5 ( 7 / + 2 7 /) + <5(7/ — 77 ) 5 ( 7 / + 77 ) = f(u). (3.12) 


3.2 Bethe Ansatz solution 


For general values of the boundary parameters k± ,£±, we have not yet succeeded to hnd 
a determinant representation analogous to ( |2.22|) of the functional relations ( |3.7| ). Never¬ 
theless, for the following two special cases of the boundary parameters, we have found such 
representations. 


3.2.1 The diagonal case k ± = 0 


The Bethe Ansatz solution for the diagonal case k± = 0 is already known !,!■ Nevertheless, 
it is instructive to see how this solution emerges from the functional relations. Indeed, when 
k± = 0, the functional relations for p > 2 can be represented as 


det 


/ A 0 
-h 
0 


V 


-hi 


' p -1 


~h'_ 1 

Ai -"•0 
—/12 a 2 


0 

0 


0 

-hi 


0 

0 


0 

0 

-K 

0 

0 


0 

0 

0 


-h. 


0 


0 

0 

0 


p —1 Ap_i 


~hn 


-ho 

0 

0 


-hi 


A r 


\ 


P~ 2 


= 0 , 


/ 


4 The last two terms of the left-hand-side were accidentally omitted in Eq. (5.2) of [h 


(3.13) 
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where 


h{u) = — sinh (u A rj) 
h'(u ) = h(—u — 2 rj ), 


2 N / , __ N sinh(2w + 2 rj) 


sinh(2u + rj) 


sinh(w + £_) sinh(« — £ + ) , 


(3.14) 

(3.15) 


and hk = h(u + rjk), h! k = h'(u A rjk), A k = A {u A rjk). We let (Qo ,Q i,... , Q p ) be the null 
vector of the matrix in (|3.13|); i.e., 


AoQo h_\Q i hg Qp 0, 

fokQk—l A A kQk hfc—lQk+l 0 j 

hp_iQo hpQp—i A A p Qp 0 . 

We make the Ansatz Q k = Q(u A rjk), where 


k = !,••• ,P~ 1 


M 


Q 0) = n sinh(u — Uj) sinh(n A Uj + rj ), 

3= 1 


(3.16) 


(3.17) 


which has the crossing symmetry Q(u) = Q(—u — rj). Eqs. ( 3.16 ) and ( [3. 15 ) imply that the 
eigenvalues are given by 


« / \ , / \ Q(« - v) , ,, \ Q( u + rj) 

A (u) = h[u) ——— -h u — rj)- 


(3.18) 


Q(u) Q(u) 

We verify that this result is consistent with both the periodicity (|3.4| ) and crossing (|3.5| ) 
properties of the transfer matrix. The requirement that A (u) be analytic at u — Uj yields 
the Bet he Ansatz equations 

h(uj) _ Q(uj + rj) 


h(—Uj - rj) Q(uj - rj) ’ 




(3.19) 


The results (p,14| ), ( |3.17| )-( |3.19|) for the transfer-matrix eigenvalues and Bethe Ansatz 
equations agree with those of Alcaraz et al. H] and Sklyanin p]. Although we have assumed 


that rj has the values ( 2.16 ), these results are true for general values of rj. As in the case of 


the closed chain, one can establish that M — A- — S z by considering the asymptotic behavior 
of t(u) for u —> oo. 


3.2.2 The nondiagonal case n + = K- = n, ^, N = odd 

Finally, we consider the nondiagonal case rt + = = n ^ 0, = ^_ = ^, N = odd. 


For this case, the functional relations again have the determinant representation (|3.13|) , with 

sinh(2u + 2 rj) 


h(u) = — sinh 2Ar (u + rj) 


sinh (2 m + rj) 


Sill 


h(u A 0 sinh(M — 0 A k 2 sinh 2 2 u) . (3.20) 
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It follows that the transfer-matrix eigenvalues and Bethe Ansatz equations are again given by 
( P.17| )-( [3.19|) . However, unlike the two cases considered earlier which have a f/(l) symmetry, 
here the value of M is fixed. Indeed, the asymptotic behavior ( |3.6| ) implies that 

M = i(JV - 1). (3.21) 

We expect that, as in the previous cases, these results hold for general values of rj. 


4 Discussion 


We have seen that an approach used by Bazhanov and Reshetikhin [IU| to solve RSOS models, 
which is based on a model’s functional relations, is also an effective means of solving vertex 
models. Indeed, we have used this method to recover the known Bethe Ansatz solutions of 
both the closed and open XXZ chains with 17(1) symmetry. Moreover, since this method 
does not rely on the existence of a pseudovacuum state, we have also been able to use this 
method to solve the special nondiagonal case ( |1.2|) of the open chain. 

Although we have focused here on vertex models associated with s/ 2 , it is clear that the 
same approach should be applicable to vertex models associated with higher-rank algebras. 

Having found a model’s functional relations, a crucial step in this method is to refor¬ 
mulate the functional relations in determinant form. We have not yet succeeded to carry 


out this step for general values of the boundary parameters of the open XXZ chain (|1 . 1| ) . 
It would clearly be useful to find necessary and sufficient conditions for the existence of a 
determinant representation of the functional relations, as well as a systematic procedure for 
its construction. We hope to be able to report on these matters in a future publication. 
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